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ABSTRACT 



This study, assuming stationary Gaussian turbulence 
model, investigates the effect on the crossing frequency of 
different spectral functions: von Karman, Kaimal and 

Teunissen for the x and z directions with five filters: 
ideal band-pass, ideal low-pass, "quadratic-type", "sine- 
type" and the Hanning. The filters have a much greater 
effect than the spectral functions. The estimated crossing 
frequency variation is as much as 50 percent among the 
quadratic-type, low-pass and sine-type filters. The Hanning 
filtering predicts crossing* rates up to thirty times, and 
the ideal one decade wide band-pass filtering predicts rates 
of eight times, higher than the ideal low-pass filtering. 
The variation, between the x and z direction, is less than 
10 percent for the von Karman spectrum, and over 40 percent 
for the Teunissen one. 
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I. INTRODUCTION 

Most flows occurring in nature and in engineering appli- 
cations are turbulent. In the earth's boundary layer, 
differential heating of the atmosphere produces pressure 
gradients, which are subsequently modified by the rotation 
of the earth, causing a complex velocity field. In this 
boundary layer the wind speed decreases as the surface is 
approached due to both the frictional drag of the surface 
and the drag of all bodies protruding into the air flow. 
These retarding forces are transmitted through the layer by 
shear forces and by the exchange of momentum due to the 
vertical movement of the air. The process of momentum 
exchange between layers is the mechanism leads to the gener- 
ation and decay of eddies which are termed turbulence. The 
resulting mixing of the air produces, along all three 
orthogonal axes, fluctuations in wind speed, commonly called 
gusts, which vary in size in both time and space. 

The effects of this atmospheric turbulence has been of 
continuing concern to the aircraft or structure designer. 
Typical turbulence related problems are: the effects of 
turbulence on the fatigue life of the structure; the 
performance of control systems in turbulence and the deter- 
mination of ultimate structural strength required to accept 
unsteady loads induced by turbulence. In an attempt to solve 
these problems, a number of statistical models of turbu- 
lence, which endeavour to describe the turbulence in terms 
of as few parameters as possible, have been proposed since 
early in this century. To analyze these problems, one needs 
to know the crossing frequency, which is the frequency with 
which a random function ( here the filtered wind speed) 
crosses a prescribed value defined by a certain averaging 
time, say the hourly mean. 
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The crossing frequency depends largely on whether or not 
the random function is filtered, the type of filter and to a 
lesser extent, the spectral density function. The effects of 
type of filter used and spectral density expression will be 
discussed in detail later. The spectral density functions 
have been proposed by Dryden [Ref. 1], von Karman [Ref. 2], 
Kaimal [Ref. 3] and Teunissen [Ref. 4], The Dryden expres- 
sion makes mathematical analysis simple but is much less 
accurate than the others. The von Karman spectral function 
was recommended by the Engineering Sciences Data Unit 
(E.S.D.U.) [Ref. 5] and this expressions was originally 
postulated for isotropic turbulence. The effect of departure 
from isotropic turbulence near the ground is allowed for by 
the variation with height and surface roughness of the 
appropriate variance O' and length scale parameter since 
they typify the intensity and size of eddies constituting 
turbulence. The Kaimal expression, which is frequently 
used, was obtained recently for the surface layer over flat, 
uniform, relatively featureless terrain in Kansas. The 
Teunissen model, which is a modified Kaimal one, was 
obtained for the generally rougher gross features of the 
upstream terrain. 

In the following discussion, we will review turbulence 
models and spectral density functions. The spectral density 
S(n), in non-dimensional form, is a function of the dimen- 
sionless parameter nL^/U, where n is frequency and L^* is the 
length scale of the turbulence and 0 is the mean hourly wind 
speed. E.S.D.U. [ Bef . 5] gives empirical expressions for 
the length scale obtained by analysing its collection of 
world-wide turbulence data. Finally, for given values of 
nL^/U, we will estimate the crossing frequencies with von 
Karman, Kaimal and Teunissen spectral expressions. 
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II. TURBULENCE MODELS 



A. DISCRETE GUSTS MODEL 

Powell and Connell [Ref. 6] defined gusts as constitu- 
ting any series of discrete velocity-time events that can be 
defined from a turbulence time series according to some 
extrinsic criterion. There are two fundamentally different 
treatments of gust time; one is that the gust time is arbi- 
trarily fixed and the other is that the gust time varies. 
Examples of this latter type of gust definition, shown in 
Figure 2.1, represent a sample of wind fluctuation in the 
form of a wind-component time series. 

The GUSTO and GUST1 models define gust events in terms 
of gust amplitudes and characteristic times. Definitions of 
amplitudes and times differ for the two models. The GUSTO 
model can be completely specified in terms of a positive 
peak amplitude + A 0 , and the time interval +T e between zero 
crossings on either side. It can also be specified in terms 
of a negative peak amplitude -A 0 and the corresponding time 
interval -T 0 . Both definitions are expected to yield 
similiar statistics. Powell and Connell combine both posi- 
tive and negative values into one set for their analysis. 

The GUST1 model can be specified in terms of a positive 
amplitude which is the peak-to-peak amplitude between 
adjacent minima and maxima, and a time interval +T, between 
the minima and the maxima; the positive sign indicates a 
positive rate of change in the variable. A comparable defi- 
nition can be made in terms of amplitude -Aj and time -T x , 
where the negative sign indicates a negative rate of change 
in the variable. 



GUST 0 Definition 




GUST 1 Definition 




Figure 2.1 Definitions of GUSTO and GUST1. 

The discrete model treats turbulence as a series of 
isolated gusts and the basic assumptions are 

1) Atmospheric turbulence can be modeled as a collection 
of isolated gusts randomly distributed along the 
structure. 

2) Gusts have random magnitude but fired shape. 

3) The structure is a deterministic linear system with 
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sufficient damping, that each encounter with one of 
these gusts results in a single significant response 
peaks. 

This model is typically used to estimate ultimate 
strength requirements, however, it is not a satisfactory one 
for all aspects of the design problem. For example, its 
application in areas such as structural fatigue, control 
system performance, or even extreme responses which involve 
lightly damped modes is highly questionable. In most 
instances turbulence does not occur as discrete gusts but as 
a continuous random disturbance. Furthermore, this model 
neglects most of the dynamic characteristics of the vehicle. 
Thus, although the necessity of evaluating the responses of 
a proposed structure to discrete gusts is still recognized 
as an important part of the design procedure, the use of the 
discrete model to calculate most response statistics has 
largely given way in recent years to the use of turbulence 
models which attempt to take both the dynamic characteris- 
tics of the structure and the continuous nature of turbu- 
lence into consideration. This is known as the continuous 
or power spectral model. 



B. CONTINUOUS NODE! 

The basic idea of continuous model is that atmospheric 
turbulence can be represented by a continuous stochastic 
process which acts as a disturbing influence on the struc- 
ture. The principal assumptions for this model are 

1} Each encounter of an structure with continuous atmos- 
pheric turbulence can be modeled as a deterministic 
linear system perturbed by three independent stati- 
onary stochastic processes, which represent the long- 
itudinal, lateral and vertical gusts components. 

2) The spectral density of each random process has a 
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form characterized by three parameters, the gust vel- 
ocity standard deviation [(S ') , the scale length (L) , and 
the averaged speed (0). 

3) Each of the three gusts components is a Gaussian 
process. 

The assumptions of the Gaussian model make it possible 
to calculate the statistics of any structure response for 
each region of turbulence as functions of the parameters L 
and O' of that region. Thus, the assumptions of structure 
linearity and the Gaussian nature of turbulence permit their 
evaluation with a minimum of difficulty. These results will 
be dependent upon the assumed values of L and o' along with 
the characteristics of the structure dynamics. 

The continuous model is used primarily to evaluate 
response statistics for selected conditions in continuous 
turbulence. For this application both the scale length and 
standard deviation of the turbulence as well as the struc- 
ture characteristics are fixed at values representative of 
certain conditions which could reasonably be expected to 
occur in service. 

The three assumptions made above need some reconsidera- 
tions. The first assumption stated that the turbulence is 
stationary and homogeneous and the dimensions of structures 
are much smaller than the scale lengths of the gusts compo- 
nents. These conditions are not always satisfied. The second 
assumption of the power spectral model is valid at low alti- 
tudes, but is not so certain at high altitudes, the problem 
made being more difficult by a lack of data at the latter 
altitudes. The third assumption stated that the turbulence 
is a Gaussian process. However, some experimental data indi- 
cate that the Gaussian turbulence model significantly under- 
estimates the frequency of occurence of high gust velocities 
[Ref. 7], Furthermore, peak gust velocity cumulative prob- 
abilities of exceedance observed in atmospheric turbulence 
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behave like exp (-x) rather than exp(-x 2 ), as predicted by 
the Gaussian distribution [fief. 8]. Recently Reeves, et al. 
[Ref. 9] have discussed in great detail a non-Gaussian 
model. On the other hand. References (6) and (10) indicate 
that filtered wind s peeds (part icularly with a band-pass 5-50 
Hz filter) is reasonably Gaussian in character. In the 
present work, the Gaussian model is used exclusively. 

The statistical quantities of most frequent interest are 
the mean, variance, probability density, autocorrelation 
function, power spectral density and crossing frequency. The 
power spectral shapes usually assumed are those proposed by 
von Karman, Kaimal and Teunissen. The crossing frequency 
with these three spectral expressions will be studied in 
more detail later. 
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III. SPECTRAL DENSITY AND PROBABILITY DEN SIT Y 
A. SPECTRAL DENSITY 

In the field of mechanical and electrical oscillations, 
a large number of phenomena are governed by equations of the 
form 



■ & 4 c * * - '«*» 

where, m, c and k are constants; f (t) is a given forcing 
function and y (t) represents the response of the system as a 
function of the time t. 

If the above time dependent force f (t) is periodic, it 
can be expressed in terms of a Fourier series; this series 
represents the sum of a large number of forces of different 
amplitudes and frequencies. All elastic bodies have natural 
frequencies, and if this force is applied to an elastic 
system, and the amplitudes of the components near a natural 
frequency are not very small, then this force will drive the 
body into resonance. The body will then go through a very 
large number of stress cycles in a short period of time, 
possibly leading to fatigue failure. 

To have stationary properties, a random signal must be 
assumed to continue over an infinite time, and in such a 
case neither the real nor the imaginary part of the Fourier 
transform converges to a steady value. In this case, we use 
the spectral density function, which has no convergence 
difficulty and which is applicable to a whole class of 
similarly generated functions. 
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Spectral functions of atmospheric turbulence provide 
information on the frequency distribution of the kinetic 
energy of the various fluctuating velocity components. Used 
in conjunction with certain transfer functions, they provide 
information about the dynamic loading on, and response of, 
building and aircraft structures in the atmospheric wind. 

A considerable number of measurements of the power 
spectra of gust velocities in the near-neutral atmosphere, 
or strong wind conditions, at varying heights and for 
different terrains are available. The von Karman spectral 
expressions, which were exclusively used by E. S.D.U. 
[Eef. 5], for each velocity component (x, y and z direction) 
are given by 



nS u (n) _ 4n u (3.1) 

<r u 2 (1 + 70.8n*) 5/6 



nS y (n) _ 4n y (1+755. 2n~y) (3.2) 

cr* { 1 + 283. 2n*)" /(S 



nSw (n) 4nn/ (1+755. 2n£) (3.3) 

c ~ ( 1+283. 2n£)"/* 



where n u = L u n/U, n v = L v n/U, n w = L^n/U and L^* is integral 
length scale which is a function only of the terrain rough- 
ness and height above ground for the x and y directions; for 
the z direction it is a function of elevation only, o'/ is 
the variance of the wind velocity fluctuations about the 
hourly mean. These equations are in common use in the engi- 
neering literature and have the advantage that the integral 
length scale are treated as ’free’ scaling parameters 

which are chosen to match the estimated scales for a partic- 
ular height and terrain type, while maintaining constant 
spectral shape. 
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Another commonly used spectral expression is the Kaimal 
one; £Ref. 3] 



nS; (n) _ 



^< 2 



0.164 ( f/f„) 
1+0. 164 (f/f„ ) S/s 



(3.4) 



where nS (n) = logarithmic power spectral density 
= variance of i 
f = reduced frequency 

f e = reduced frequency at the intercept of the 
extrapolated inertial subrange slope with the nS 4 * (n) / 
<$l =1 line; or in our notation, f©; =0.041z/L t *. 

The Kaimal expression is mere common in the meteorolo- 
gical literature and was obtained as a best-fit to surface- 
layer measurements over uniform, flat, relatively 
featureless terrain. 

Teunissen [Ref. 4] observed that von Karman model- 
appears to be somewhat better than the Kaimal model at low 
frequencies. The Kaimal model largely underestimates the 
spectral content at these frequencies for all compo- 
nents (x, y, and z) . On the other hand, the Kaimal model shape 
appears to be better than that cf the von Karman model. 

Teunissen has proposed modified Kaimal spectral expres- 
sions, given by 



nSti (n) _ ' 1 0 5f 

u, 2 (0.44 + 33 

nSy (n) _ 1 7f 

u^ 2 (0.3 8+9.5f/ /3 

n5„(n) _ 2 f 

u* 2 0.4 4+5.3f 5/a 



(3.5) 



(3.6) 



(3.7) 
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where z=nL^/U and u* is friction velocity. These expres- 
sions were derived from the data which were obtained from a 
horizontal array of tower-mounted propeller anemome- 
ters (z=11m) during a five-hour period for which the mean 
wind direction was virtually perpendicular to the main span 
of the array, in the city of Toronto, Canada. 

As for the reason for the -higher spectral energies found 
by Teunissen in the low-f reguency region of the spectra in 
comparison with the Kaimal values which was measured in 
Kansas, this is not entirely obvious, apart from the rougher 
terrain for Teunissen experiment. Teunissen has suggested 
that the power spectra over a wide range of terrain types 
may possibly be represented by using above expressions with 
appropriate values of such ’terrain scaling’ parameters. 

B. PROBABILITY DENSITY 

In general, a knowledge of the spectral density of a 
random process does not enable us to determine its prob- 
ability density or distribution. However, a common assump- 
tion which is reasonable in many applications is that 
atmospheric turbulence is a ’normal* or Gaussian process, as 
discussed above, with a probability density function for 
which 



ta) =Trir exf, [ _ii/2ir * 2 ] 



(3.8) 



where i=u(t), v(t) or w{t) are wind-speed fluctuations, 
which have a zero mean. Note that Eguation (3. 8) depends only 
on the standard deviation of the gust velocity. Thus, if the 
probability distribution of o' is known, the distribution 
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function of the gust velocity can be found by integration. 



1 • 0 • f 



P = f p (i) di (3.9) 

00 

In practice, atmospheric turbulence contains patches of 
a significantly non-Gaussian nature (particularly in the 
lower 30m) when larger gusts and longer lulls may occur more 
frequently then indicated by the Gaussian distribution. 
Kaimal et al. £Ref. 10] observed that the assumption of a 
stationary Gaussian process has validity only for the 
filtered data, not for the unfiltered time series. 
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IV. CROSSING FREQUENCIES 



The instantaneous horizontal wind speed is considered to 
be composed of a mean longitudinal component 0 plus fluctu- 
ating component v^* , where the subscript i designates the x, 
y, or z component in a coordinate system with x oriented 
along the mean wind vector and z upwards. Then the wind 
speed as a function of time is given by 



The average frequency of positive-slope level crossing 
Noi , with which the wind speed fluctuations exceed their 
zero value, for a Gaussian distribution of turbulence, is 
given by 



where n c is cut-off frequency, i is u, v,or w and F (n) is a 
filter function. 

The value of N 0x - depends on the mean wind speed U , L f * and 
turbulence intensity <%/U , which is a measure of the magni- 
tude of turbulence fluctuations and defined as the ratio of 
the standard deviation of the instantaneous fluctuating 
velocity component to the mean wind speed averaged over 
arbitrary time, usually hourly, and depends only on the 
height z and the terrain roughness z e , in a stable 
atmosphere. 



V (t) = U + v;{t) 



( 4 . 1 ) 




( 4 . 2 ) 
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A. FILTER FUNCTIONS 



It is clear that the quality of the crossing frequency 
estimates obtained by means of Equation (4. 2) will depend 
strongly on the form of the filter function F(n) and, hence, 
on the selection of the time function h(t). Thus we will 
consider some examples of the function h(t) and filter func- 
tion F{n) corresponding to them. Let the time function be 
given by 



h L (t) = 2 sin (27Ct/T) (4.3) 

T 2TCt/T 

If we take Fourier transform, the Equation (4. 3) becomes 



F x (n) =/■ 1 
1/2 
l 0 



when | n | < 1/T 
when |n| =1/T 
when |n| >1/T. 



(4.4) 



Graphs of the functions h^ (t)' and F^ (n) are shown in Figure 
4.1 and Figure 4.2 respectively and they represent an ideal 
band-pass filter. 

The simple impulse response function h 2 (t) is the ordi- 
nary moving average; 



h 2 (t) =r 1/T when -1/2 <t< +T/2 (4.5) 

1 0 when | t | > T/2 

The corresponding filter function F 2 (n) (F 2 (n)=H 2 (n) 2 ), where 
H 2 (n) is the Fourier transform of the function h 2 (t) , has 
the form 



F 2 (n) 



sirntTn 

TCTn 



(4.6) 
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TIME FUNCTION 




Figure 4.1 Time Function h^ (t) , h 2 (t) , h 3 (t) and h^.(t) 



The function h 2 (t) and F z (n) are shown in Figure 4. 1 and 
Figure 4.2 respectively. The figure shows that F 2 (n) is a 
low-f reguency filter with cut-off frequency n c =1/T. The 
shape of this filter shows that frequencies above n c are 
transmitted nonuniformly and those below n c are completely 
suppressed. 

The impulse response function h 3 (t) , called the smoothed 
moving average is: 



hj (t) =M+cos (2TC/T) t when 
T 

0 when 



-T/2 <t< T/2 
1 1 1 > T/2 



(4.7) 
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FILTER FUNCTION 



in 




Figure 4.2 Filter Function Fj (n) , F 2 (n) , F 3 (n) and F^ (n) . 



The corresponding transfer function F 3 (n) has the form 



f 3 (n) 


= f sirncTn^ 


Tt 4 ' 




(4.8) 




V TCTn ) 


[it 2 - {TTTn) 2 ] Z 




The advantage 


of F 3 (n) 


over F 2 (n) , 


is 


that it has much 


smaller side maxima. 








Another commonly used 


filter function 


is the guadratic 



filter function given by 

( n ) = 1 (4.9) 
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One can find the impulse function of Equation (4 . 9) by 
inverse Fourier transform, i.e.. 




(4. 10) 



The second term of above equation is zero, since, it is an 
odd function. Thus Equation (4 . 10) becomes 



where K 0 is the modified Bessel function. 

The functions h^ (t) and (n) are shown in Figure 4.1 
and Figure 4.2. 

The cut-off frequencies are defined formally as the 
frequencies at which the filters reach the first zero. But, 
when a real filter is used, we must find some equivalent 
cut-off frequency, which is defined as the cut-off frequency 
of an equivalent ideal filter whose area is equal to that of 
the real filter. [Bef. 11] Greenway [Ref. 12] and Powell 
and Connell [Bef. 6] have used the value at the half- 
power (3db) point of the filter as the cut-off frequency. 
Vinnichenko et al. [Bef. 11] defined the equivalent cut-off 
frequency as 




or 



h^(t) =^K 0 (n c t) 



( 4 . 11 ) 



n, 



c 




o 



(4. 12) 
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EFFECTS OF CUT-OFF FREQ. 




Figure 4.3 Effects of Cut-off Frequencies. 

According to the Equation (4. 1 2) , the cut-off frequencies 
depend on the filtering function. The effects of cut-off 
frequency on the crossing frequencies are shown in Figure 
4.3, which shows that the higher the cut-off frequency leads 
to lower zero crossing rates. 
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B. CALCULATION OF CBOSSING FREQUENCIES 



In this section the frequency of zero crossings are 
estimated using von Karman, Kaimal and Teunissen spectral 
density functions with various filter functions which we 
have discussed above. 

1 . Osina the von Karma n Spect ral De ns i ty F unct ion 

Introducing the von Karman spectral density function 
for the x-direction S u (n) into Equation (4. 2) yields 



N 



ox 



r ~ nc * 

P n J F (n) dn 


/ f nca F(n) dn 


J(1+70.8nJ) 5/ V 

> HC| 


J (1+70.8n*f 4 
n cl J 



1 Vz 



For numerical calculation, set y=n/n C| and D=n C2 /10n C! 
Equation (4. 13) transformed to 



(4. 13) 



and 



~ n ci 



ploD 

J.O + 



y 2 dy 



iop 



70 . 8n y 3 f 




ay 

70. 85* 



-i'/z 



(4. 14) 



Y 3 ) 



S/6 



where N 0 * is the dimensionless frequency of positive- si ope 
level crossings for the x direction N ox L u /U. Equation (4. 14) 
is for the ideal band-pass filter function of D decades 
filter width. 

The upper integral of Eg. (4.13) will not converge as 
n c2 goes to infinity because the empirical expression for 
the turbulence spectrum employed does not account for the 
high-freguency range where viscous dissipation damps out 
turbulence fluctuations. One can achieve convergence by use 
of a suitable filter function. Introducing Equation (4 . 9) 
into Equation (4. 13) yields 



N 



ox 



n c 


00 

f y <3y 


/r d y 1 




^ { 1+ 70 . 8 n c 2 y 3 ) 5/ ^ ( 1 ♦ y 2 ) y 


' J(1 + 70.8n 3 y J f*(1 + y 3 ) 



(4.15) 
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If we take F 2 (n) , here called the sine filter function, 
instead of quadratic filter, with the cut-off frequency 
n =0.44/1 which Greenway [Ref. 12] has used, and substitute 
into Equation (4. 6) , we find 

F 2 (n)=sin 2 -(0.44nn/n c )/(0.44it n/n c f. (4. 16) 

Introducing Eq ua tion (4. 1 6) into Equation (4. 13) and setting 
y=0.44un/n c yields 



N oa = n c 



0. 44TC 




sin z (y) dy 



J/l 



70. 8 (n c /0.44X) 2 y 2 ] 
[sin 2 (y) /y 2 ] dy 
[l+70.8(n c /0.44X) 2 y 2 ] 






(4. 17) 



% 



If we use Fj (n) , which is the Hanning frequency 
window, and equate cut-off frequency n c =0.375/T as derived 
by Vinnichenko et al. [Ref. 11], then substitution into 
Equation (4. 8) gives 



F 3 (n) = 



sin (0 . 3757Cn)\ It 

0. 375]tn /[id 2 - (0.375irnf] Z 



(4. 18) 



Substituting Equation (4 . 1 8) intc Equation (4. 13) and setting 
y=0.375nn/n c and =n c /0. 37 5TL yields 



r »° .2 
I sin 

J 4(1 + 70.8*1 * 



(Y) dy 

y z f /6 






y 2 ) 



£ 



[sin 2 (y) /y 2 ] dy 



(1+70.8^/ y 2 f/* (U 2 -y 2 ) 2 



(4. 19) 
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TABLE I 

Crossing Frequencies (von Karman x- direct ion) 



lie 


IDBP 


QUAD 


SINE 


HANNING 


0. 100 


0.3972 


0. 1215 


0.0865 


17. 4480 


0. 125 


0.4838 


0. 11 16 


0.0955 


21.8074 


0. 150 


0.5709 


0. 1250 


0. 1 1 17 


26. 1670 


0. 175 


0.6585 


0. 1381 


0. 1232 


30.5267 


0.220 


0.8172 


0. 1604 


0. 1424 


38.3741 


0. 500 


1. 8193 


0. 2737 


0.2391 


87.2049 


1.000 


3. 6235 


0.4298 


0. 3699 


174. 4062 


2.000 


7. 2394 


0.6765 


0.5744 


348. 8093 


5. 000 


18.093 1 


1. 2365 


1.0383 


872.0234 


10.000 


36.1847 


1.9555 


1.6267 


1744.0481 


15.000 


54.2766 


2. 5585 


2. 1230 


2616.0669 


22.00 0 


79.6055 


3.2989 


2.7307 


3836.9021 


50.000 


180.9212 


5. 6928 


4. 6950 


8720.2461 


75.000 


271. 38 16 


7.4555 


6.1414 


— 


100.000 


36 1. 8425 


9. 0288 


7.4326 


— 


150.000 


54 2. 7637 


11.8268 


9.7292 


— 


175.000 


633.2241 


1 3 . 1 0 54 


10.7787 


— 


200.000 


723.6853 


14.32 43 


11.7794 


— 


225.000 


8 14. 1465 


11.4933 


12.9392 


— 


250.000 


904.6067 


16.61 96 


13.6641 






nc =cut-off 


frequency 








IDBP=ideal 


band- pass 


filter 






Q0AD= ’’quadratic-type 


" filter 






SINE= ’’sine- 


-type’’ fii 


ter 






HANNI NG= HANNING filter 






IDLP=ideal 


low-pass 


filter 





The results of numerical calculations of Eqs.(4.14), (4.15), 

(4.17) and (4.19) are given in table I. The plots for above 
equations are shown in Figure 4.4. The plots showing how 
well these crossing frequencies co-incide depend on the 
length scale, mean wind speed and filter functions. 

The von Karman spectral density function for the 
z-direction is given by Equation (3. 3) . The equation of the 
crossing frequency for the ideal band-pass filter function 
can be obtained by substituting Equation (3. 3) into 
Equation (4. 2) and is given by 
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CROSSING FREQUENCY(N 0 *) 



Vi 



N02 



=n 



ci 




|oI> 

(1+755. 2n C i y 2 ) y z ay 
(1+283. 2 ng y 2 )” /6 
?1+755.2n ci 2 y a )dy 
(1+283. 2n c , 2 y 2 )" /6 



(4. 20) 



VON KARMAN(X-DIRECTION) 




Figure 4.4 Crossing Frequencies with von Karman 
Spectral Expression for the x-direction. 
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Using the same procedure as for the x-direction, the 
z-direction crossing frequency for the quadratic filter 
function is found to be 



N 



oz 





°°( 1+755. 2n/y 2 ) y 2 dy 
(1+283. 25^* f* (l+y 2 ) 
*°( 1+755. 2n c a y 2 ) dy 
( 1+283. 2n^y 2 )^(1+y a ) 



l /2 



(4.21) 



and, for the sine filter function, if one defines 
n c /0. 44TC. =*[ c the crossing frequency for the z-direction 
becomes 



N 



OH 




f [1*755.2 (q c y f J sin 2 (y) dy 
(1 + 283. 2 n?Y 2 )" /6 
r Il+755. 2 (>f c y) 2 J[sin 2 (y) /y^dy 
J 0 (1 + 283. 2\Iy 2 ) , ‘ / (> 




(4. 22) 



Introducing filter function F 3 (n) into Equation (4. 20) and 
setting y=0.375nn/n c and n c /0. 375it= r ( c yields 



f f [l +755. 2 (*f fc y ) 2 ] sin 2 (y) dy 
N ez =\ J (1 + 283. 2^cV) % (Tl a -y 2 ) 2 






‘ [1+755. 2 (yf/y^JC sin 2 (y) /y*]dy 
(1+283. 2^ c 2 y 2 ) //< (U Z -y a ) 2 



(4.23) 
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TABLE II 






Crossing Freguencies (von Karman z-direction) 



Ik 


IDBP 


QUAD 


SINE 


HANNING 


0. 100 


0.38 12 


0. 121 


0.0816 


17. 4448 


0. 125 


0.4686 


0. 136 


0.0934 


21.8044 


0. 150 


0.5567 


0. 150 


0. 1044 


26. 1641 


0. 175 


0.6454 


0.163 


0.1144 


30.5238 


0. 220 


0. 8059 


0. 185 


0. 1320 


38.3716 


0. 500 


1.8137 


0. 297 


0.2188 


87.2036 


1. 000 


3.6207 


0. 452 


0.3375 


174.4053 


2.000 


7. 2380 


0.696 


0.5238 


348. 8093 


5.000 


18.0925 


1 .248 


0.9458 


872. 0229 


10.000 


36. 1 844 


1. 951 


1.4870 


1744.0461 


15.000 


54.2764 


2. 535 


1. 9407 


2616. 0693 


22.000 


79.6054 


3. 246 


2.4978 


3836.8989 


50.000 


180.9214 


5.493 


4.2994 


8720.2266 


75.000 


271. 3823 


7. 127 


5.6256 


— 


100.000 


361. 84 30 


8.627 


6. 8094 


— 


150.000 


542.7634 


11.129 


8.9150 


— 


175.000 


633. 2231 


12. 330 


9. 8773 


— 


200.000 


723.6838 


13.474 


10.7948 


— 


225.000 


814. 1440 


14.572 


1 1. 6748 


— — 


250.000 


904. 6057 


15.63C 


12. 5228 


— 



The results of numerical evaluations of Egs. (4.20), 
(4.21) ,(4.22) and (4.23) are listed in the Table II. The 
plots for above equations are shown in Figure 4.5. 

For comparison, the crossing frequencies with von 
Karman spectral expression for the x and z directions are 
plotted in Figure 4.6. The crossing frequency for the z 
direction is less than the x direction except ideal band- 
pass filter case, in which the frequencies are identical for 
both directions. 

The crossing frequency plots presented in Figures 
4.4 - 4.6 offer a detailed view of how the crossing 
frequency depends on the filter function. Thus, in designing 
a system, the type of the filter has to be chosen carefully 
to match the behavior of the particular system. 



32 



CROSSING FREQUENCY(N oz ) 



VON I< A RM AN(Z-DI RECTION) 




Figure 4.5 Crossing Frequencies with von Karman 
Spectral Expression far the z-direction. 



2. U sing the Kaimal Spectral Density Funct ion 

From Equation (3. 4) , the asymptotic spectral behavior 
can be written as 

nS<- (n) _ f (f/f 0 ) V3 ; f»f 0 (4.24a) 

<^ 2 1 0.164 (f/f e ) ; f<<f 0 (4.24b) 
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VON KARMAN(X AND Z DIRECTION) 




Figure 4. 6 Crossing Frequencies with von Karman 
Spectral Expression for the x and z-direction. 



The expression for the low frequency behavior of the 
spectrum is given 

= ( 4 U/t)f 

(4.25) 

where I t * represents the length scale. Equating 

Equation (4 . 25) and Equation (4. 24b) becomes 

L<- = 0.041z/ (f.)* 



(4. 26) 



For the x-direction only 
L u = 0. 04 Iz/fo-x 
or f 0 y = 0. 04 1z/L u 

Substituting Equation (4. 27) into (3.4) yields 



(4. 27) 



nS u (n) _ 0. 164 (f/Q. 041z/L u ) 



or 



Ou 

nS u (n) 



1+0. 164 (f/0. 04 1z/L u ) 
4n 



u 



1+33. 641n 



^3 



(4. 28) 



where n = nL u /U. Then the crossing frequency for the 
X-direction becomes 



No* = 



f- 

in 



F (n) n dn 



(1 + 33.641n 5/ 3) / J(1+33.641n 5/ 3) 



f- 

J{ 1 + 






F (n) dn 



(4.29) 



Substituting the quadratic filter 
Equation (4 .29) yields 

.oo „ Vz 



function into 



A/ 

N 



OX 



I- 

-4(1+ 



n 2 dn 



33. 641n g/ 3) (n* + n z ) 



f- 

J(U 



dn 



J (1 + 33.641n^3) (n 2 + n 2 ) 



(4. 30) 



Replacing n/n c by y in Equation (4.30) leads to 

1 Vz 



N o* =n c J 



r 

-/[i* 



y dy 



33. 641 (n^ y)^3] (H-y^) 



£ 



dy 



[l + 33. 641 (n c y) r/3 ](1 + y a ) 



(4.31) 
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in Equation (4. 31) gives 



Setting y=e 



No* =n c 




oo 



e 3 * dx 



[1+33.641 (npe^)^] (1+e zX ) 



Vz 



/, 



e* dx 






(1+33.64 1 (n c e x ) 5/ 3j (He 2 *) . 



If we introduce the sine filter function 
Equation (4 . 29) , the crossing frequency becomes 



N 



OX 




sin 2 (0 . 44rcn/n c ) n 2 dn 
(1+33.641n ff /3) (0. 44it/n c ) 2 
00 sin 2 (0 . 44^/^ ) dn 

(1+33. 641 n ff /3) (0 . 4 4Tt/n c ) 2 



,Vl 




Substituting 0.44jtn/n c for y in Equation (4. 33) , and 
n c /0 . 44Tt=l'| c , yields 



N 



ox 



=n c 




sin 2 (y) dy 
1+33.641 (He y ) 5/3 
" ’[sin 2 (y) /y 3 ] dy 
1 + 33.641 (5f c y) 5/3 



V 2 



For the Hanning filter, substituting Equation (4 . 1 8) 
Equation (4. 28) , setting y=0 .375mn/n c and 

n c /0 . 375H=>( C , yields 



(4.32) 



into 



(4.33) 



define 



(4.34) 



into 

define 
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(4.35) 









f- 

J.U* 



sin 2- (y ) dy 
[1+33.641 (^ c Y) y/3 3 



l / 2 



(n; 2 -y 2 ) 2 



u 



Csin 2 (y) /y 2 J dy 



[1+33. 641 (?[ c y) 7/ 3] (it 2 -y 2 ) 2 



TABLE HI 

Crossing Frequencies (Kaimal x-direction) 



Sc 


IDBP 


£240 


SINE 


HANNING 


0. 100 


0.4056 


0. 134 


0. 0868 


17.4532 


0. 125 


0.4938 


0. 152 


0. 1001 


21.8126 


0. 150 


0.5819 


0. 169 


0.1126 


26. 1720 


0. 175 


0.6702 


0. 184 


0. 1243 


30.5316 


0. 220 


0. 8295 


0.210 


0. 1441 


38. 3790 


0.500 


1.8305 


0. 340 


0.2445 


87.2088 


1.000 


3. 6321 


0.518 


0.3811 


174.4094 


2.000 


7. 2456 


0.799 


0.5943 


348. 8120 


5. 000 


18.0968 


1.434 


1.0744 


872.0224 


10.000 


36.1872 


2. 249 


1.6879 


1744.0486 


15.000 


54.2786 


2.931 


2.1987 


26 16.0667 


22.000 


79.6070 


3. 769 


2.8249 


3836. 9033 


50.000 


180.9221 


6. 471 


4.8291 


8720. 2461 


75.000 


271.3823 


8.464 


6.2923 


— 


100.000 


36 1. 8423 


10.243 


7.5860 


— 


150.000 


542.7634 


13.324 


9.8619 


— 


175.000 


633. 2239 


14.733 


10.8962 


— 


200.000 


723.6846 


16.102 


11. 9076 


— 


225.000 


814. 1450 


17.414 


12.8333 


— 


250.000 


904.6052 


1 8.678 


13.7017 


— 



The results of numerical calculations of Egs. (4.29), 
(4. 32), (4.34) and (4.35) are given in Table III and the 

plots for above equations are shown in Figure 4.7. 

It is clear from the Eg. (3.4) and Eg. (4.26) that, 
for the Kaimal spectral density function, the difference in 
crossing freguencies between the x and z directions is only 
in the length scale . Thus the plots for the z-direction 
are exactly same as the x-direction. 
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Figure 4.7 Crossing Freguenci.es with Kaimal 
Spectral Expression far the x- direct ion. 

3. Dsinq the Teunissen Spectral Density Fun ctio n 

The modified Kaimal spectral density expressions, 
which were derived by Teunissen, are given in Egs. (3.5), 
(3.6) and (3.7) . These eguations were derived by adding a 
scaling parameter to each of the Kaimal spectral eguations 
and evaluating these parameters from the observed data. 
These eguations are function of friction velocity u^. . In 
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order to estimate the crossing frequency as we have done 
earlier, we need to change that equations as a function of 
variance. In E.S.D.U. [Ref. 5], the relations between 0^ 



and u^ are defined as 

0.4O2 /u*=(<£/U) log (z/Zo) (4.36) 

Teunissen derived the average values of q2 / u^., given by 

<T w /u* = 2. 84 (4.37a) 

o^/u* = 1.92 (4.37b) 

<52,/u* = 1.27 (4.37c) 

From Equation (4. 36) and Equation (4. 37) , one can get 

uj (u-component) = 0.124 (4.3 8a) 

u£ (v-component) = 0.271 (4.38b) 

u* (w-component) = 0. 620 < 5 "^ (4.38c) 



and substituting Eg. (4.38) into Eqs.(3.5), (3.6) and (3.7) 

yields 



nS u (n) 



1 3n 



(0. 44 + 33n) 5 ^ 3 



(4.39) 



nSy (n) 
^ 2 



4. 61 n 



(0.38 + 9.5n) 5/3 



(4.40) 



nS w (n) _ 1 . 24n 






0. 44 + 5 . 3n^ 3 



(4.41) 



Using methods similar to the above, the crossing frequency 
for the x-direction, with the ideal band pass filter, 
becomes 
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< 0 * 



r n ca n 2 dn j r nca dn 
J (0 .44+33n?7 J ( 0. 4 4 + 33 n) 5/a 

^ n c» Hci 



*/2 



or 



N 



ox n ci 



r " c> y^dy /r^L 

J (O.V>»33n c , yf 3 J(.0.‘ 

«ci ' nci 



dy 






ci yf 3 /J(0.44 + 33n cl yf /3 
nci 



(4. 42) 



If one uses the sine filter function, with cut-off 
frequency n c =0.44/T, the crossing frequency becomes 



N 0 * = 



f- 

Jo (0 



sin* (0. 44un/n c ) n 2 dn 
/ Q (0. 44+33n) 5/>3 (0 .44fln/n c ) 2 



I 



sin 2 (0.4 4nn/n c ) dn 
(0.44+33n) 5/3 (0.44nn/n c ) 3 



V? 



(4.43) 



Setting y=0.44itn/n c into Eg. (4.43), and define n c /0 . 44Tl=i , ( c , 
yields 






N o x = n c 



r 00 

r 


sin 2 (y) dy 


1 r[ sin 2 (y) /y 2 ] dy | 




(0.44+3 3^ c y) % j 


/I (0.44 + 33^ c y)H 


; filter function. 


N 05( becomes 


1 


r °° ? 2 dy 




l/ i 


J, 


j (0. 44+3 3fT c y)^ 


(Uy 2 ) 






r dy 




L J 


I (0.44+33n c y) S/ * 


—4 

+ 





l /2 



(4. 44) 



(4.45) 
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TABLE IV 

Crossing Frequencies (Teunissen x-direction) 



Sc 


I DBF 


£UAD 


SINE 


IDLP 


0. 100 


0.3722 


0.060 


0.0425 


0.0334 


0. 125 


0.4623 


0.069 


0. 0494 


0.0394 


0. 150 


0.5534 


0.078 


0.0559 


0.0449 


0. 175 


0.6439 


0.086 


0. 0620 


0.0502 


0.220 


0.8068 


0.099 


0.0724 


0.0590 


0. 500 


1.8201 


0.169 


0. 1251 


0 . 1044 


1.000 


3.6294 


0.266 


0. 1983 


0. 1672 


2. 000 


7.2478 


0.420 


0.3142 


0.2664 


5. 000 


18. 1031 


0.768 


0. 5754 


0.4915 


10.000 


36. 1952 


1.213 


0.9083 


0.7801 


15.000 


54. 2873 


1.581 


1. 1846 


1.0219 


22.000 


79.6162 


2.039 


1.522 8 


1.3188 


50.000 


180.9321 


3.521 


2. 6002 


2.2786 


75.000 


271.3931 


4.612 


3.3801 


2.9851 


100.000 


361. 8535 


5. 536 


4. 0689 


3.6157 


150.000 


542.7749 


7. 266 


5.2816 


4.7372 


175.000 


633.2354 


8.053 


5. 8554 


5.2496 


200.000 


723.6 956 


8.804 


6. 2870 


5. 7380 


225.000 


814. 1565 


9.525 


6.8032 


6. 2065 


250.000 


904.6 169 


10. 219 


7.301 1 


6.6579 



The results of numerical calculations of Eqs. (4.42) , 
(4.44) and (4.45) are given in Table IV and the plots for 
these equations are shown in Figure 4.8. 

The modified Kaimal spectral expression for the 
z-direction is given by Eg. (4.4 1). If we put ideal band pass 
filter function and spectral expression for the z-directon 
into Eq. (4.2) it yields 
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CROSSING FREQUENCY(No*) 



or 






N 



O £ 



= n, 



y\ci 



ci 



k 

.nc, 



y2dy 



44+5.3 (n cl yf /3 




r\c2 



dy 



2 

(4. 46) 



0. 44+5. 3 (n cl y) 

n C i 



% 



TEUNISSEN(X-DIRECTION) 




Figure 4.8 Crossing Frequencies using Teunissen 
Spectral Expression for the x-direction. 
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If one introduces the sine filter function, the equation for 
the crossing frequency becomes 



N 



05 



Ij 

f- 

1 (c 



sin 2 (0 . 4 4itn/nc ) n 2 dn 
(0 .44+5. 3n g / a ) (0. 4 4ttn/n c ) 2 

sin 2 (0. 44itn/nc) dn 

(0.44+5.3n ?/ 3) (0. 4 4un/n c ) 2 / 






or 




He 




2 

sin (y ) dy 
0. 44+5. 3 (n c Y) y/3 
[sin 2 (y) /y 3 j dy 
0.44+5. 3 (»[ c y) 5/3 




If one introduces the quadratic filter function, the 
tion for the crossing frequency becomes 



/-.t* 



N oz " n c< 



J o [0.44+5.3 (n c y) 5/3 j (1 + 



Y 2 ) 



if 



dy 



44+5. 3 (n c Y) S/i ] (1 + Y 2 ) • 



Vi 



The results of numerical calculations of Eqs. ( 
(4.48) and (4.49) are given in Table V. The plots for 
equations are shown in Figure 4.9. 



4.47) 



4. 48) 



eg ua- 



(4.49) 



.46) , 
above 
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TABLE V 

Crossing Frequencies (Teunissen z-direction) 



Sc 


IDBP 


£UAD 


SINE 


IDLP 


0. 100 


0.4454 


0.178 


0. 1098 


0.0559 


0. 125 


0.5367 


0.201 


0.1258 


0.0691 


0. 150 


0.6263 


0.222 


0. 1408 


0.0818 


0. 175 


0.7151 


0.241 


0.1551 


0.0941 


0.220 


0.8740 


0.274 


0.1794 


0.1153 


0. 500 


1.865 1 


0.438 


0.3047 


0. 2268 


1.000 


3.6560 


0.660 


0.4758 


0.3783 


2.000 


7.2610 


1. 007 


0. 741 4 


0.6104 


5.000 


18. 1053 


1.789 


1.3353 


1 . 1232 


10.000 


36. 1926 


2.791 


2.0917 


1.7725 


15.000 


54. 2827 


3.629 


2.7242 


2.3144 


22.000 


79.6 102 


4.658 


3.5000 


2.9785 


50.000 


180.9 239 


7.973 


5. 9955 


5. 1216 


75.000 


271 .3835 


10. 4 10 


7.8197 


6.6982 


100.000 


361. 8437 


12. 591 


9.4420 


8. 1050 


150.000 


542.7644 


16.470 


12.3095 


10.6064 


175.000 


633. 2249 


18. 242 


13.6158 


11.7493 


200.000 


723.6853 


19.9 32 


14.8606 


12.8388 


225.000 


814. 1458 


2 1 . 5 53 


16.0474 


13.8837 


250.0 00 


904.6062 


23. 1 15 


17. 1816 


14. 8905 



For comparison, the crossing frequencies with the 
von Karman, Kaimal and Teunissen spectral expressions for 
the x-direction with various filter functions are plotted in 
Figure 4.10. As shown in Figure 4.10, the crossing 
frequencies predicted by the Kaimal spectral expression are 
higher than for the others. In Figure 4.11, the crossing 
frequencies, with all three spectral expressions and the 
z-direction, are compared. The figure characteristics are 
similar to those for the x-direction, except that the 
crossing frequencies are lower. 
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Figure 4.9 Crossing Freguencies with Teunissen 
Spectral Density Function for the z-direction. 



C. EXAMPLE CALCULATIONS 



In the previous section of this study, we derived equa- 
tions for the crossing frequency. The purpose of this 
section is to illustrate the application of these equations 
to a given turbulence condition, and to estimate the 
crossing freguencies for a low -pass filter. 
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VON KARMAN, Kf IMAL& TEUNISSEN(X) 




Figure 4.10 Comparison of von Karman. Kaimal and Teunissen 
Spectral Density Function for the x-direction. 

The dimensionless crossing frequency, N 0 ^ , is a function 
of n c (=nl^ /D) . The length scale parameter and surface 
roughness parameter z c are given by E.S.D.O. [Ref. 5]. For 
convenience these data are included here as Figure 4.12 for 
length scale parameters, and Figure 4.13 for the surface 
roughness parameters. 
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Figure 4.11 Comparison of von Karman, Kaimal and Teunissen 
Spectral Density Function for the z-direction. 



Suppose we have a site for which 
z = 20m 
n t = 50Hz 

0 = 25m/sec and the turbulence condition is that 
for an off-sea wind in a coastal area. 

With this choice of turbulence condition, the surface 
roughness parameter z 0 can be found in Figure 4. 13, and is 
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z o =0.00 1. The length scale parameters are given in Figure 
4.12, and are given as L u =105, l v =36 and L w = 6.5. 

The next step is calculation of the dimensionless 
cut-off frequencies, i.e., 

n cu = nL u /0 =210 (4. 50) 

n cv = nL v /U = 72 

n ca/ = nL w /U = 13 

The crossing frequencies are calculated by the use of 
the given crossing frequency vs cut-off frequency graphs, 
with dimensionless cut-off frequencies which are given in 
Eg. (4.50) . 

The final step is to convert the dimensionless crossing 
frequency to a dimensional frequency by multiplying the 
factor U/L^ , i.e., 

N 0 * = N 0 * 0/I u (4.51) 

Noe = N oZ U/l w 

1 The results of calculation cf the crossing frequency for 
a given turbulence condition are given in Table VI. 
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TABLE VI 








Example 


Calculations 


of the 


Crossing Frequency 






von Karman ( x-d irec 


tion) 






REMARK 


IIEP 


ID IP 


quad 


SINE 




Nox 


75.0 


10.2 


14.6 


12.0 




Nox 


17.857 


2. 428 


3.476 


2. 857 






{S'cu =210. 


0, U/Lu=0. 238) 








von Karman ( z-direction) 






REMARK 


I EBP 


ID IP 


QUAD 


SINE 




Noz 


4.70 


1. 50 


2.40 


1. 80 




Noz 


18.076 


5. 769 


9.231 


6.923 






( n cw — 1 -3 / 


U/Lw=3 


.846) 








Kai mal (x and z -direction) 






REMARK 


I EBP 


I DIP 


QUAD 


SINE 




Nox 


75.0 


10.9 


16.50 


12.2 




Nox 


17.857 


2. 595 


3.928 


2. 905 




Noz 


4.70 


1. 70 


2.65 


2.05 




Noz 


18.076 


6. 538 


10. 192 


7.885 






Teunissen (x 


and z- 


direction) 






REMARK 


I EBP 


I DIP 


22 AD 


SINE 




Nox 


75.0 


6. 08 


9.30 


6.50 




Nox 


17.857 


1. 448 


2.214 


1. 548 




Noz 


4.70 


2. 10 


• 3.30 


2. 40 




Noz 


18.076 


8. 076 


12. 692 


9.230 


. 
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Figure 4.12 Values of length Scale Parameter L; 
(Reproduced from EsDO 74031). 
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Figure 4.13 



Values of the Surface Roughness Parameter z„ 
(Reproduced from ESDO 740J1). 
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V. CONCLUSIONS 



This study has estimated the crossing frequency based on 
the three spectral density functions, which assumes turbu- 
lence to be a homogeneous, stationary Gaussian process, with 
five different filter functions for the x and z directions. 

This estimation has shown, for the x direction, that the 
Kaimal spectral function gives about one percent higher 
crossing frequency than the von Karman, and forty percent 
higher than the Teunissen one, at ^ = 100 for quadratic-type 
filtering; for the ideal low-pass and sine-type filters, the 
difference between the Kaimal and von Karman spectral func- 
tion is only a few percent, tut for the Teunissen one is 
more than forty percent than for the other spectrums. 

For the z direction, the Teunissen spectral function 
gives about ten percent higher crossing frequency than the 
Kaimal one and about twenty percent higher than the von 
Karman one, at n c =100 for the quadratic-type, sine-type and 
low-pass filtering. However, this difference becomes to zero 
as non-dimensional cut-off frequency goes to zero. 

The crossing frequencies for band-pass and Hanning 
filtering are only a fraction of a percent different for the 
three spectral functions. 

The present study shows that Hanning and band-pass 
filtering predicts higher crossing frequencies than low-pass 
filtering. This difference is greatly amplified with the 
higher frequency. 

If some site statistics are available, as for the 
example calculation, then the number of crossings that the 
wind-speed fluctuations exceed their zero value per a 
certain period can be computed and an estimation of fatigue 
life can be obtained. 
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As a consequence of the results studied herein, it is 
suggested that, in future work, careful consideration should 
be given to the type of filtering used to model structural 
behavior and determine more accurately the fatigue life of a 
system. 
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